Mass Dependence of Ultracold Three-Body Collision Rates 
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We show that many aspects of ultracold three-body collisions can be controlled by choosing the 
mass ratio between the collision partners. In the ultracold regime, the scattering length dependence 
of the three-body rates can be substantially modified from the equal mass results. We demonstrate 
that the only non-trivial mass dependence is due solely to Efimov physics. We have determined 
the mass dependence of the three-body collision rates for all heteronuclear systems relevant for 
two-component atomic gases with resonant s-wave interspecies interactions, which includes only 
three-body systems with two identical bosons or two identical fermions. 
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The achievement of quantum degeneracy in ultracold 
gases with different atomic species |l| has driven exper- 
imental studies of several novel phenomena. The ob- 
servation of interspecies Feshbach resonances in boson- 
fermion mixtures j^] allows considerable flexibility for ex- 
ploring new regimes by controlling the interspecies inter- 
actions. Boson- mediated Cooper pairing 3], for instance, 
can substantially increase the critical temperature for a 
phase transition to the Bardeen-Cooper-Schrieffer (BCS) 
regime. The collapse of the fermionic component as 
well as the phase separation of both components can 
be studied along with the creation of ultracold po- 
lar molecules 0. It is worth noting that heteronuclear 
boson-fermion molecules are composite fermions which 
allows a new type of crossover between an atomic Bose- 
Einstcin condensate and a molecular Fermi-type super- 
fluidity (6|. Heteronuclear boson-boson and fermion- 
fermion mixtures have also been studied , but no Fesh- 
bach resonances have yet been reported. 

The magnetic field sensitivity of the hyperfine states 
is the key to controlling the interatomic interactions in 
ultracold gases. By applying an external magnetic field 
near a diatomic Feshbach resonance, the s-wave scatter- 
ing length a, which characterizes the low-energy inter- 
atomic interactions, can take any value from the weakly 
(a — > 0) to the strongly (|a| — ^ oo) interacting limits. 
Even though two-body loss processes can usually be min- 
imized by using resonances in the lowest hyperfine states, 
three-body loss processes can still be substantial. Fortu- 
nately, near the resonance, when \a\ ^ rg (with tq being 
the characteristic range of the interatomic interactions) 
processes such as vibrational relaxation, X+X2 ^ 
three-body recombination, X + X -\- X X + X2 , and 
collision-induced dissociation, X + X2^X + X + X, no 
longer depend on the details of the interactions and uni- 
versal predictions can be made. 

Recent experiments have underscored the importance 
of knowing the a dependence of three-body rates in or- 
der to determine the atomic and molecular lifetimes. In 
fact, three-body losses have been used to locate Fesh- 
bach resonances 01 and to create ultracold molecules 
. While general results for threshold and scattering 



length scaling laws of three-body equal mass systems [loj 
have been obtained, however, there are no similarly gen- 
eral scaling laws for heteronuclear systems. The specific 
case of recombination in a two-component Fermi gas has 
been investigated, though, and found to scale as for 
a > 0, and minima were predicted as a function of the 
mass ratio between the collision partners • 

In this Letter we demonstrate that the mass ratio has 
a large impact on ultracold three-body coUisional losses, 
allowing a certain degree of control. Using the simple 
physical picture developed in Ref. 0], extended to in- 
clude heteronuclear systems, we have determined that 
the scattering length scalin g la ws can differ substantially 
from the equal mass results . For instance, in a system 
with two identical fermions that are much heavier than 
the third atom, relaxation of weakly bound heteronu- 
clear molecules scales approximately as — an even 
stronger suppression than the a~^'^^ scaling found when 
all three atoms have equal mass 0| . This scaling was 
derived in Ref. to explain the long molecular lifetimes 
observed experimentally for molecules formed of fermions 
in different spin states. It was their long lifetimes that 
made further experiments with these molecules feasible 

lllj ] . The stronger suppression found here might open 
other experimental avenues. 

In our picture, the mass dependence enters via the 
three-body effective potentials and is associated with Efi- 
mov physics 01 . In fact, we demonstrate here that the 
only non-trivial mass dependence is due solely to Efimov 
physics. We determine the mass dependence of the col- 
lision rates for three-body systems relevant to all two- 
component atomic gases with resonant interspecies s- 
wave interactions, which includes systems with two iden- 
tical bosons or two identical fermions. 

For short-range two-body interactions, the ultracold 
behavior of three-body systems can be derived from 
three-body effective potentials and couplings ^3 which, 
in the adiabatic hyperspherical representation, are deter- 
mined from the adiabatic equation [isj . 

H,4R, r!)$,(i?; Q) = f/,(i?)$,(i?; n). (1) 

This equation is obtained from the full Schrodinger equa- 
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tion by fixing the hyperradius R, leaving only dynam- 
ics in the hyperangles ft through /fad- By expanding 
the total wave function on the adiabatic basis the 
Schrodinger equation (in atomic units) is reduced to: 



1 

2^1 dR^ 



E 



Vuu'Fy' = EFy 



(2) 



where E is the total energy, Fi, is the hyperradial wave 
function, and v \s a, collective index that represents 
all quantum numbers necessary to label each channel. 
In the present case, the three-body reduced mass is 
pL—mj ^(5((5 -I- 2) where to is the mass of the distinguish- 
able particle and b is the ratio of the different atom mass 
to the identical atom mass. This equation describes the 
collective radial motion under the influence of the effec- 
tive potential Wy as well as any inelastic transitions con- 
trolled by the nonadiabatic couplings V^yi . 

In the limit \a\ ^ tq the effecti ve p otentials Wy depend 
crucially on Efimov physics [HIlJ. In Ref. Il0|, we de- 
scribed a scheme that uses this fact to classify all equal- 
mass three-body systems with resonant s-wave pairwise 
interactions. In this scheme, each system falls into one 
of two categories: those with an attractive dipole (i?~^) 
potential in the range ^ i? <^ \a\ and those only with 
repulsive dipole potentials. Which category a particu- 
lar system falls into depends on its symmetry J'^ (total 
angular momentum J and parity tt) and the identical 
particle permutational symmetry, i.e., whether they are 
bosons or fermions. For heteronuclear systems, this ba- 
sic classification scheme still holds, and the extension of 
the analysis of the three-body rates requires only slight 
modifications which, in turn, come almost entirely from 
the 8 dependence of the dipole potential strengths. 

For each category, there are three distinct regions in R 
that characterize the effective potentials (for a schematic 
picture, see Fig. 1 in Ref. lCr|): R < ro, ro<Ci?^|a5|, 
and R->\af\ where ai^\yj 8{8 + 2)/((5 + l)]i/^a. Replac- 
ing a by in these definitions is the first mass-dependent 
modification of our previous analysis. In the asymptotic 
region, _R3>|a5|, the potentials can be derived analyti- 
cally pjij. They are associated with molecular chan- 
nels, which represent atom-molecule scattering, and with 
three-body continuum channels, which represent colli- 
sions of three free atoms, and are given respectively by 



= E^v^ \ J and = ^ ^ ' — —. (3) 
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2/^i?2 



The molecular bound state energy E^ii is labeled by the 
ro-vibrational quantum numbers v and V] I is the atom- 
molecule relative angular momentum; and A labels the 
eigenstates of the hyperangular kinetic energy. 

For ro<Ci?<C|a5|, the potentials for the molecular and 
continuum channels [Eq. ^] are modified due to Efimov 
physics, establishing our classification scheme. In the 



first category, an attractive dipole potential occurs in the 
highest vibrationally excited s-wave molecular channel 
for a > and in the lowest continuum channel for a < 0. 
The potentials for all higher-lying channels are repulsive. 
These potentials are conveniently parametrized by the 
coefficients sq and s^: 



Wy{R) 



2^i?2 



and Wy{R) 



"v 4 
2^i?2 



(4) 



In the second category, the potentials for the weakly 
bound molecular channel and the continuum channels are 
repulsive and parametrized by coefficients po and p^: 



Wy{R) 



and Wy{R) 



2^iR^ 



(5) 



In both categories, however, deeply bound molecular 
channels are essentially independent of a. In the above 
equations, the coefficients sq, s^, po, and Pi, depend on 
the number of resonant pairs, the number of identical 
particles, and the mass ratio S between collision partners 
|l4| . For i? < To, the potentials depend on the details 
of the interatomic interactions and can lead to resonance 
effects due to three-body shape or Feshbach resonances. 

In the adiabatic hyperspherical representation, inelas- 
tic transitions at ultracold temperatures proceed via 
tunneling in the initial collision channel to where the 
coupling Vyyi peaks — the coupling, of course, drives 
the transition. We have previously shown that a sim- 
ple WKB approximation to the tunneling probability 
through the potentials given by Eqs. ©-© is sufficient 
to determine both the threshold and the scattering length 
scaling laws for equal mass three-body collisions JOJ . For 
heteronuclear systems, as, sqj s^, pq, and generate the 
only non-trivial mass dependence in the three-body rates 
and follow directly from Efimov's analysis. 

For systems and symmetries from our first category 
[with an attractive potential Q], relaxation for a>0 is 



sinh(27/) 



sin~[so ln(a5/ro) + <I>]+sinh (77) 



and for a<0, Vrcioc^'-i£;^^iiro'+^ Ecoii=E-E„i' is the 
collision energy, $ is an unknown short-range phase |lfll | , 
and 77 labels parameters related to the inelastic transition 
probability at small distances . These parameters can 
depend on the masses nontrivially and can lead to reso- 
nance effects, depending upon details of the interactions. 
Recombination for a>0 and a<0 is given by 



-,21+1 



; (6) 
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K^OQu^^^E^ — ^ — , \as 

sin^ [so ln(|a5|/ro) + $]+sinh-'(?7) 



2A+4 
S ' 

2A+4 

(7) 
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For systems and symmetries from our second category 
[with a repulsive potential ISJ], relaxation for a > is 



(8) 



while it is V^oi /i' ^-E^oU^o'^^ f'^'' a < 0. Recombination 
for a > and a < is given by 



as as 
K3 0,^i'~'E\ro/\as\)'P°\as\'^+\ 



,2A+4 



(9) 



The mass dependence of sq and po can substantially 
modify the a dependence of the three-body rates. The 
coefficients and p^, however, do not affect the scaling 
laws. Figure ^ shows sq and po for systems with two 
identical bosons {BBX, 6—mx/mB) and with two iden- 
tical fermions {FFX, 5—mx/mp). These coefficients are 
determined analytically following Ref. after proper 
symmetrization. It is important to notice that for J > 
the effective potentials can change from attractive to re- 
pulsive and vice-versa, changing the category in which 
the system is classified, and thus modifying its collisional 
properties. It happens, for instance, for 2+ BBX sys- 
tems at 6c « 0.0259 [Fig. [Ha)] and for 1" FFX systems 
at (5c « 0.0735 [Fig. [Tib)] . 
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FIG. 1: Mass dependence of so and po, for (a) two identical 
boson and (b) two identical fermion systems. 

For bosonic systems BBX , the dominant contribution 
for relaxation and recombination is J"^ = 0+ with I = 
and A = 0, and is described by Eqs. ®-tZI) due to the 
presence of an attractive dipole potential. In this case, 
the variations of Sq with S [Fig.^a)] changes the number 
of Efimov states. These variations manifest themselves 
in the three-body rates through the locations of the min- 
ima and peaks in Eqs. ©-(0). We note that the mass 
dependence in this case does not modify the scaling laws. 

In contrast, for fermionic systems FFX, the power law 
behavior does change with mass for the dominant contri- 
butions to relaxation ( = 0+ with / = 0) and recombi- 
nation {J^ = 1~ with A = I). These changes follow from 
the mass dependence of the repulsive dipole potential due 
to Po [Fig. nib)] combined with Eqs. Q-®. It follows 
that relaxation for FX + F collisions scales as a^^'^P" for 



a > and is suppressed for all 6 since 2<po<4 [Fig.Q^b)], 
approaching its greatest suppression, a"*", as (5 and 
its least, a~^, as (5 ^ cxo. Suppression of FX + F colli- 
sions can thus be much stronger than a~^'^^, found for 
S — 1 FF' + F collisions of fermionic atoms in differ- 
ent spin states ^1] ■ Extremely long-lived heteronu- 
clear molecules in ultracold boson-fermion and fermion- 
fermion mixtures might thus be possible. The modifica- 
tions to the recombination scaling law occurs for a < 
and S > 0.0735 where it scales as a^~'^P° with < po < 2 
[Fig.^b)], otherwise it scales as a^, leading to an asym- 
metry between a > and a < 0. The strongest asym- 
metry is obtained when (5 — > oo, where recombination 
scales as and |ap, respectively, so that greater colli- 
sional stability against recombination might be expected 
for a < 0. 

As mentioned above, recombination for FFX systems 
with S > 0.0735 was studied in Ref. 11] and found to 
scale as a^ for a > which agrees with our prediction. It 
was also predicted that the recombination rate oscillates 
as a function of S with zeros at S ~ 0.0735, 0.1160, and 
as ^ ^ oo (note that the mass ratio in [l^ is 1 /6) . These 
zeros were ascribed to a decoup ling of the three-body 
continuum and bound channels [ijj, leading to elastic 
scattering only, and have also been interpreted as inter- 
ference effects (17j . In our formulation, however, no such 
zeros are predicted. In fact, in this range of i5, the poten- 
tials for both continuum and bound channels, Eq. Q, are 
repulsive and interference effects are expected to be sup- 
pressed by tunneling. Further, we have noticed through 
representative numerical calculations that the continuum 
and bound channels never decouple and that the nona- 
diabatic couplin g p eaks at i? « rg and i? « a, in accord 
with our model 'lO|. While we have observed minima in 
our calculations, they do not correspond to the minima 
discussed in Ref. ^3] ■ They are most likely a consequence 
of using a two-body potential with finite a, and we expect 
them to disappear in the limit a — > oo considered in [111 ] . 
Preliminary numerical results support this expectation. 

TableQlsummarizes the threshold and scattering length 
scaling laws for the three-body rates, including dissoci- 
ation D3. The three dominant partial waves (as deter- 
mined by their energy dependence) are shown for each 
process. The table thus indicates the terms expected to 
be important for finite energies or \a\ 00. For BBX 
systems, these contributions are not expected to be im- 
portant for relaxation with a > or recombination with 
a < 0. For FFX systems, however, the higher partial 
waves are comparatively more important. 

In two-component ultracold atomic gases with reso- 
nant interspecies interactions, the important three-body 
processes are those that involve both atomic species, as- 
suming that intraspecies interactions are not resonant. 
In this case, there are only two relevant systems and the 
competition between the collision rates for each system as 
well as the density of each species determines the atomic 
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TABLE I: Threshold and scattering length scaling laws for 
three-body rates. For 2+ BBX and 1~ FFX systems the 
rates are given for 5 < 5c and 5 > Sc. Boldface indicates the 
leading contribution at threshold. 



process, we have demonstrated that the mass dependence 
in three-body collisions is intimately related with Efimov 
physics. In bosonic systems, the mass dependence af- 
fects the features due to Efimov physics but leaves the 
dominant scattering length scaling laws unchanged. In 
fermionic systems, however, the mass dependence modi- 
fies the scattering length scaling laws substantially from 
the equal mass results. The stronger suppression found 
for relaxation of weakly bound heteronuclear molecules in 
boson-fermion mixtures, suggests the possibility of long- 
lived heteronuclear fermionic molecules. 

This work was supported by the National Science 
Foundation. 



and molecular lifetimes. In boson-fermion mixtures, for 
instance, relaxation for BF + F collisions decreases with 
a while for BF + B it increases. In this case, if all bosons 
are bound in BF molecules, the molecules are expected 
to be long-lived. Otherwise, the molecules will be rapidly 
quenched by collisions with the bosonic atoms. 

Vibrational relaxation, as we have shown, can be con- 
trolled by choosing the atomic species according to their 
mass and consequently the molecular lifetime can be sub- 
stantially modified. We show in Table ITll the suppression 
predicted for BF + F collisions of commonly used alkali 
atoms. The table includes systems used in recent ex- 
eriments: ^SNa-^Li and ^■^Rb-'*°K 0, as well as "^Li-^^Li 
We have included H due to the prospects for using 
it for sympathetic cooling Molecule-molecule colli- 
sions also contribute to the molecular lifetime but 
in boson-fermion mixtures the molecules are composite 
fermions so that only p-wave molecule-molecule collisions 
occur. In this case, molecule-molecule collisions are sup- 
pressed for ultracold temperatures, and atom-molecule 
collisions are expected to be dominant. 
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TABLE II: Scattering length scaling laws for relaxation (a > 
0) oi BF-\-F collisions in boson-fermion mixtures based solely 
on mass. 

In this Letter, we have explored the mass dependence 
of the ultracold three-body coUision rates for systems 
with two identical particles that interact resonantly with 
a third particle. We have derived the energy and scat- 
tering length dependence for all three-body inelastic col- 
lision rates that apply to two-component ultracold quan- 
tum gases with resonant interspecies interactions. In the 
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